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Abstract 

In this paper, we study the relationship between the network-based inference method 
and global ranking method in personal recommendation. By some theoretical analysis, 
we prove that the recommendation result under the global ranking method is the limit 
of applying network-based inference method with infinity times. 
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1. Introduction 


Personal recommendation [1] has been a hot topic recently because of its wide and 
effective application in business, and many recommend methods have been developed. 
For example, global ranking method (GRM), network-based inference (NBI) method 
i3i (also called the mass-diffusion (MD) method or probabilistic spreading (ProbS) 


algorithm), heat-spreading (HeatS) algorithm 


(also called the heat conduction 


(HC) or heat diffusion process). Many new algorithms have been proposed based on 
these methods, see |5|- 7| and references therein. 

In this paper, we will concentrate on the relationship between NBI and GRM in the 
theoretical view. By studying the NBI, we propose a new algorithm by using the NBI 
method with multiple times. Using the matrix analysis technique, we rigorously prove 
that the result on GRM equals to that under the limit of NBI with infinity times. 
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2. GRM and NBI: A review 


The following definitions and statement of GRM and NBI have been given in 
here we rewrite it to keep the self-integrity of this paper. 


A, 


2.1. GRM 

At first, we state the recommendation process of GRM. 

GRM sorts all the objects in the descending order of degree, and recommends those 
with highest degrees. 

GRM lacks of personalization, but it is widely used since it is simple, the well 
known “Yahoo Top 100 MTVs”, “Amazon List of Top Sellers”, as well as the board of 
most downloaded articles in many scientific journals, can be all considered as results 
of GRM. 


2 . 2 . NBI 

NBI sets up a bipartite to discriminate the object-set and the user-set, then uses the 
diffusive idea to determine the importance of node i in node j’s sense. 

At first, we present the definition of bipartite. 

Definition 1. (See jSJ) A simple graph G is called bipartite if its vertex set V can be 
partitioned into two disjoint sets V\ and V2 such that every edge in the graph connects 
a vertex in V\ and a vertex in V2 (so that no edge in G connects either two vertices 
in V\ or two vertices in V2). When this condition holds, we call the pair {V\, V2) be a 
bipartition of the vertex set V of G. 

Consider a general bipartite network G(Vj, V2, E), V\ = {01, ■ • • , o„ } means the 
object-set, and V2 = {u\, 112, ■ ■ ■ , u m } means the user-set. The nxm adjacent matrix 
A is defined by A = (an), where an = 1 if user it; has already collected object Oi and 
an = 0 otherwise. Moreover, assume the initial resource located in the z-th object is 
f(oi ) > 0 . A sketch map of the bipartite network with three objects and four users is 
given in Figure [T] 

The first key point for NBI is that: the resource in an arbitrary V\ node should be 
equally distributed to its neighbors in V2, here we use notation Vj —> V2 to denote this 
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Figure 1: Sketch map of the bipartite network with initial resources x,y,z 


step. After this step, all the resource in V\ flows to V2, and the resource located on the 
^-th V2 node is: 


f(ui) = 

where k(oj) is the degree of object j. 
result under step V\ —> V2. 


sr' a jif(°j) m 

h k ^ 

Figure U gives a sketch map of the diffusive 
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x/3+y/2+z/3 y/2+z/3 


x/3+z/3 


Figure 2: Sketch map of the diffusive result under step V\ —> V 2 
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The second key point for NBI is that: the resource in any V2 node should be equally 
distributed to its V\ neighbors, here we use notation V2 —» V\ to denote this step. After 
this step, all the resource in V2 flows back to V\, and the final resource located on o, is 


f'iPi) = E 
1 = 1 


aufjui ) 
k(ui) 


m n p / \ 

\' an ^^ajij(Oj) 

Kui) ^ k( 0j ) 


n 


3 =1 


where 


( 2 ) 


Wij 


1 

k{Oj) 


E 


i =1 


Q'ilQ'jl 

k{ui) 


( 3 ) 


Therefore, if we denote F^> = ,f'(o n )) T , W = ( w^j ) £ f? nxra anc | 

F = (/(oi),-" ,/(o„)) T , then 


F (t) = wjr 


( 4 ) 


For the example in FigureQ]and Figured after this step, see Figured the resources 
at these three nodes are denoted by x', y', z', which can be calculated as 
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Figure 3: Sketch map of the diffusive result under step V 2 —> Vi 
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NBI recommendation algorithm: The recommendation algorithm aims at pre¬ 
dicting user ui s personal opinions on those objects ui has not yet collected, l = 
1 , • • • , m. Set the initial resource located on each node Oj, j = 1 , • • • ,n of Vi as 

f(oj) = aji (6) 

That is to say, if the object Oj has been collected by u/, then its initial resource is unit, 
otherwise it is zero. The initial resource can be understood as giving a unit recom¬ 
mending capacity to each collected object. Therefore, 

n n 

f(°i) = H w n a n ( ? ) 

j=i j=i 

42 For any user m, all his uncollected objects o :] are sorted in the descending order of 

43 and those objects with highest value of final resource are recommended. 


44 3. Some theoretical preparations 


45 At the first glance, one cannot find any relationship between GRM and NBI, be- 

46 cause GRM does not consider the personal difference, while NBI considers it; more- 

47 over, numerical examples show that NBI is better than GRM in personal recommenda- 

48 tion | 2 j. In order to bridge the gap between them from the theoretical point, we should 

49 first present some useful lemmas and properties. 


Property 1 . For the matrix W = ( Wij ), where Wij is defined in 0, it is column 
normalized, that is to say 

n 

£«* = 1, j = 1, • ,n- (8) 

i= 1 

Proof: According to the definition of wij , in order to prove the above property, we 
just need to prove that the following conclusion holds: 




§ k (°j) IrJ 


= 1. 


Vj = 1, • 


( 9 ) 
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Because k(oj) is independent on parameters i and Z, so we let k{of) = p(j), where 
p(j) is an integer depending on j, which means that there are p(j) nodes in V 2 con¬ 
necting j. Without loss of generality, we assume their indexes are Zi, I 2 , ■ ■ • , 
Therefore, with the definition of adjacency matrix, the left of © can be represented as: 




Q'il&jl 


7~[ k (°f) T^i k ^ 


EE 


Qil&jl 


Qi,l. 


( a i,li , a i,l2 _|_, 

p(j)~i\ k (uh) K u h) k (ui pW ) 


E a i,h , ST' 


0*1,12 


E 


Oil 


pU) 


pU)[^K u i J ~i k ( u h) ~i k ( u iru))- 

= —— (1 + 1 ”1 -= 1 

PW ---' 

pU) 


50 The proof is completed. 


51 Remark 1. In fact, it is easy to obtain the above property from the point of matrix. At 

52 first, from the definition ofWij, we have 

W = {AU- 1 ){0~ 1 A) T , (10) 

where W = ( u>ij) £ R nxn , diagonal matrix U = diag{fc(rti), • • • , k(u m )} £ R mxm 
means the degree matrix in V 2 , so AU^ 1 is a matrix with its each column sum is 1, i.e., 

(I) - " 1 l) 1X n ■ {AU~ X ) = (1, ■ • • ,l) lxm , 

Similarly, the diagonal matrix O = diag{/c(oi), • • • ,fc(o„)} £ R nxn means the de¬ 
gree matrix in set Vj, so O^ 1 A is a matrix with its each row column is 1, i.e., (O -1 A) T 
is also a matrix with its each column sum is 1. Therefore, 

(l,-- - ,l) lxro -(0- 1 A) T = (l,-. - ,l) lxn . 

53 Thus, 

(1, • • • , l)lxn • W = (1, • • • , l)lxn • (AU-')(0-'A) T = (1, • ■ ■ , l)lxn, 

54 That is to say, matrix W is also a matrix with its each column sum is 1. Furthermore, 

55 one can also get that matrix W has a left eigenvector (1, • • • , 1) T £ R nxl correspond¬ 
ed ing to eigenvalue 1. 
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Moreover, in the form of matrix, we can also easily get one right eigenvector for 
matrix W is (fc(oi), • • • , k(o n )) T £ R nxl corresponding to eigenvalue 1, since 


W ■ (Ar(oi), — ,k(o n )) T 

=(AU~ 1 ){0~ 1 A) T ■ (fc(oi), • • • , k{o n )) T = ( AU~ 1 )A t O~ 1 • (fc(oi), • • • , k{o n )) T 
={AU~ 1 )A T • (1, • • • , l)£ xl = AU~ 1 (k(ui), ■ • ■ , fc( Mm )) T 
=-4(1, ■ ■ • . l)™ xl = (fc(oi), • • • , k(o n )) T 


For example, in Figure[I] the adjacency matrix is A = 


110 1 
0 110 
0 111 


W = {AU- 1 ){0~ 1 A) T 

( 


= A ■ 


1 0 
0 1/3 
0 0 
0 0 


0 

0 

1/2 

0 


A 2 


^ 1 1/3 0 1/2 ^ 

0 1/3 1/2 0 

^ 0 1/3 1/2 1/2 y 

f 11/18 1/6 5/18 ^ 

1/9 5/12 5/18 

^ 5/18 5/12 4/9 j 



V 0 

0 

1 1/3 

0 

o N 

1/3 

1/2 

1/3 

0 

1/2 

1/3 

V 1/3 

0 

1/3 ) 


, thus 


1/3 0 

0 1/2 


0 
0 

1/3 ) 


( 11 ) 


57 Therefore, we obtain the matrix W, which can also be found in (0), all the elements 

58 are nonnegative, and its each column sum is 1. Moreover, (1,1,1) T is the left eigen- 

59 vector for W corresponding to the eigenvalue 1, while (3, 2,3) T is the right eigenvector 

6 0 for W corresponding to the eigenvalue 1 . 

6 1 Next, we will discuss some properties of this matrix W. In order to do that, some 

62 useful lemmas should be introduced. 
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Lemma 1. (See Jsj/) (Gersgorin Disc Theorem) Let A = () £ R nxn and 

n 

G' j {A)= 

denote the deleted absolute column sums of A. Then all eigenvalues of A are located 
in the union ofn discs 

n 

IJ{~ e C : \ z ~ a n\ G Cj}. 

i =i 

Lemma 2. (See (Perron-Frobenius Theorem) For n x n matrix A, which is irre¬ 
ducible nonnegative matrix, then 

1. p(A) > 0, where p(A ) is the spectral radius; 

2. p(A) is an eigenvalue of A; 

3. there is an vector x > 0 and Ax — p{A)x; 

4. p(A) is an algebraically (and hence geometrically) simple eigenvalue of A. 

Assumption 1. In the following discussion, we will always assume graph G is con¬ 
nected. 

In fact, this assumption is not very strong. Because if G is not connected, then 
G can be split into at least two connected components, for example G i and G 2 , then 
V\ = Vj 0 Vf and V 2 = V.] 0 l/ 2 2 , such that (V^ 1 , V^ 1 ) and (Vj 2 , V^ 2 ) are two bipar- 
tites with no connection, and this case contradicts to our aim of using resource location 
feedback. Moreover, under this condition, one can also know that matrix W is irre¬ 
ducible. 

4. Main results 

Now, it is time to bridge the gap between GRM and NBI. Recalling the process 
of NBI, one can find that the two steps V\ —» V 2 and V 2 —> V± can be regarded as a 
round of resource diffusion, and the obtained resource vector is determined by the 
original resource vector F and the matrix W, see ©. Inspired by NBI, we propose a 
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multiple rounds of resource diffusion algorithm, which can be described as 


Vi -► V 2 -► Vi ->• V 2 Vi ->• V 2 -> Vi ->■ V 2 -A Vi -)■ • • • 

| 1st round | 2nd round | 3rd round | 4th round | (12) 

p ^(!) ir(2) jr(3) f (4) ... 

After IV round, the resource location ft'W j s determined as: 

jfW = WF (Ar_1) = W 2 F < ' N ~ 2 ' 1 = ■■■ = W n ~ 1 F ( ' 1 '> = W N F. (13) 

In the next, we will explore the property of matrix W and W N , in order to finalize 
the relationship between NBI and GRM. 

Theorem 1. For the matrix W € R nxn defined in ( 1701 ), suppose it it irreducible, then 
it has the following properties: 

1. M must have an eigenvalue \\ = 1, and its multiplicity is 1. Suppose ei and e r 
are the corresponding left eigenvector and right eigenvector, and eje r = 1. 

2 . e r > 0, and e r can be chosen as: 

e r = (fc(oi), k(o 2 ), ■■■ , k(p n )) T (14) 

while 

n 

ej = a(l, • • • , l) T ,a = l/^fc(oj). (15) 

j=i 

3. The other n — 1 eigenvalues Xj £ C, j = 2, 3, ■ ■ • ,n satisfying \Xj \ < 1. 

4. lim W N = e r eT. 

iV-H-oo 1 

Proof: According to the Gersgorin disc theorem (Lemma|T}, Perron-Frobenius theo¬ 
rem (Lemma|2]i, Property 1 and Remark 1, one can easily get the conclusions 1, 2 and 
3. Next, we will concentrate on proving the fourth point. 

Denote J as the Jordan form of W, i.e., 

W = PJP~ \ (16) 


9 


93 

94 

95 

96 

97 

98 

99 
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where J can be written in the form as: J = diag{l, J 2 , ■ • • , J q }, where 

( X 1 \ 

A 1 


Jb = 


V 


,0 = 2 


A 1 

X J 

where A are chosen from A», i = 2, ■ • • , n, so \Re(Xi)\ < 1, and 

->■ 0, IV —>• + 00 . 

Therefore, if we denote IF* = lim IF :V , then 

N —>- + oo 

W N = (PJP" 1 )^ = PJ N P~ l Pdiag(l, 0, • • ■ , 0)P" 1 , N + 00 , 
i.e., W* = Pdiag(l, 0, • • • , 0)P _1 . 

Because AP = PJ, so the first column of P is e r - Similarly, because P~ X A = 
JP _1 , so the first row of P -1 is ej. Since P -1 P = /, therefore, e[e r = 1, which is 
satisfied in the first point. 

In all, we can get that W* = e r ej. The proof is completed. 

For example, for the matrix W in ( 111 b . simple calculations can show that eigen¬ 
values of W are: Ai = 1 ,A 2 = 0.4034, A 3 = 0.0689; while e r = (3,2, 3) T and 
ei = 1/8 • ( 1 , 1 , l ) 2 for eigenvalue 1 , so: 


1 3 N 


OO 

OO 

00 

OO 

CO 

OO 

2 

•1/8-(1,1,1) = 

2/8 2/8 2/8 

V 3 J 


CO 

OO 

CO 

OO 

CO 

OO 


W* = e r e{ = 


Next, we use the Matlab to verify the correctness of our claim. 


(17) 



1 0.4691 0.2870 0.3395 N 


1 0.4129 0.3392 0.3609 ^ 

w 2 = 

0.1914 0.3079 0.2701 

;W 73 = 

0.2262 0.2727 0.2587 


v 0.3395 0.4051 0.3904 y 


^ 0.3609 0.3881 0.3804 } 


1 0.3903 0.3606 0.3693 N 


( 0.3812 0.3692 0.3727 ^ 

w 4 = 

0.2404 0.2591 0.2536 

;W 5 = 

0.2461 0.2537 0.2514 


^ 0.3693 0.3803 0.3771 y 


v 0.3727 0.3772 0.3758 j 
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104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 



1 0.3775 0.3727 0.3741 ^ 


( 0.3760 0.3741 0.3746 ^ 

w 6 = 

0.2484 0.2515 0.2506 

; W 7 = 

0.2494 0.2506 0.2502 


v 0.3741 0.3759 0.3753 y 


y 0.3746 0.3754 0.3751 y 


( 0.3754 0.3746 0.3749 ^ 


^ 0.3752 0.3748 0.3749 ^ 

w 8 = 

0.2497 0.2502 0.2501 

;W 9 = 

0.2499 0.2501 0.2500 


v 0.3749 0.3751 0.3751 y 


y 0.3749 0.3751 0.3750 y 


( 0.3751 0.3749 0.3750 ^ 


( 0.3750 0.3750 0.3750 ^ 

w 10 = 

0.2500 0.2500 0.2500 

; W 11 = 

0.2500 0.2500 0.2500 


v 0.3750 0.3750 0.3750 y 


^ 0.3750 0.3750 0.3750 y 


Obviously, W N can converge to W* in (fTTb when TV = 11, therefore, our claim is 
correct. 


Based on the above theorem, considering the limit case W*, for any user /, since 
f'iPi) = £" =1 w* 0 f{oj) = i w ij a ji, therefore, we have 

( /'(oi) \ ( ciii \ [ an \ ( k(oi) \ 


f(o 2 ) 




= w* 


«2 1 


^ Uni J 


T 

= e r e l ■ 


« 2 1 


^ Uni J 


= ak(yi)e r = ctfc(t/;) 


k{o 2 ) 

\ k (°n) ) 

(18) 


i.e., for recommendation one only needs to see the value of e r , while e r denotes the 
vector composed of degree, and this is just the GRM. 


Remark 2. The use of left eigenvector and right eigenvector in recommendation sys¬ 


tems can be retrieved to / j/. Moreover, it is also widely adopted in the analysis of 


synchronization and consensus literature, see 




Remark 3. In fact, using the Hamilton-Cayley Theorem, matrix W satisfies the poly¬ 
nomial /(A) = |A/— with degree n, i.e., W n +aiW n ~ 1 + •• -+a n -iW+a n I = 0. 
Moreover, according to the so-called minimal polynomial for matrix W, there exists 
a polynomial with degree n' < n, such that P(W) = W n + a\W n _1 + ■ • • + 
a n '-\W + a n 'I = 0. Therefore, for any integer N, from the theory of polynomial 
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division, there exist a quotient Q{W) and a remainder R{W), such that W N = 
Q(W)P(W) + R(W), where the degree of the remainder is less than the degree of 
the divisor P(W), that is to say, any W N ,N = 1, 2, • • ■ can be equally described by a 
polynomial with degree less than n'. In [7], the authors consider the second and even 
the third round, and simulations show that they exhibit better recommendation result 
than just only one round. Since we just want to emphasis the relationship between NBI 
and GRM from the theoretical view, here we omit the numerical simulations. Interested 
readers are encouraged to investigate this problem. 

5. Conclusion 

In this paper, we first extend the NBI which can be regarded the first round of 
resource diffusion to the multiple rounds of resource diffusion algorithm. Then by 
rigorous theoretical analysis, we finally prove that the GRM is just the limit of our pro¬ 
posed algorithm. That is to say, for the multiple rounds of resource diffusion algorithm, 
NBI is the case with N = 1, where N means the times of diffusion, while GRM is the 
case with N = +oo. We bridge the gap between GRM and NBI successfully. 
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